Bosenova and three-body loss in a 85 Rb Bose-Einstein condensate 
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Collapsing Bose-Einstein condensates are rich and complex quantum systems for which quanti- 
tative explanation by simple models has proved elusive. We present new experimental data on the 
collapse of high density 85 Rb condensates with attractive interactions and find quantitative agree- 
ment with the predictions of the Gross-Pitaevskii equation. The collapse data and measurements 
of the decay of atoms from our condensates allow us to put new limits on the value of the 85 Rb 
three-body loss coefficient K% at small positive and negative scattering lengths. 
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While most experiments with dilute gas Bose-Einstein 
condensates have employed atomic species with repul- 
sive interactions, it has long been known that interesting 
and exotic physics manifests in attracting condensates. 
These include macroscopic quantum tunnelling [1 , the 
formation of soliton trains and vortex rings [U [3J, and 
the violent collapse and explosion known as the 'bosen- 
ova' [3J-[B] . The first evidence for the collapse of attract- 
ing Bose-Einstein condensates was found by Sackett and 
coworkers, who analysed the thermal equilibration of a 
sample of 7 Li atoms with negative scattering length that 
was cooled below the critical temperature [7 . Soon af- 
ter this work, condensate collapse was directly observed 
in pioneering experiments at JILA which revealed a 
host of interesting dynamics and prompted a surge of 
theoretical interest [SWTS']. More recently, the collapse of 
dipolar chromium BECs has been observed, displaying 
the striking d-w&ve symmetry of long-range dipole-dipole 
interactions in excellent agreement with theory [SJ. 

However, while initial mean-field analysis of the JILA 
experiment using the Gross-Pitaevskii (GP) equation was 
able to qualitatively account for most of the experimental 
observations, including the formation of atomic 'bursts' 
and 'jets' [8 , 12 -15], further investigation exposed a quan- 
titative discrepancy between theory and experiment in 
the time taken for the condensates to collapse [TT] . This 
was especially puzzling as the short time, low density 
phase of the experiment is exactly where the GP equa- 
tion should be an excellent approximation. This disagree- 
ment, of about 100%, could not be eliminated by more 
complex quantum field calculations [T51 ITFJ] . and has led 
to the development of competing models for the collapse 
mechanism [17_. Yet amid the extensive theoretical work 
on this phenomenon that has continued in recent years, 
there has been a notable absence of further experimental 
data, and the discrepancy between theory and the 85 Rb 
experiment remains unresolved. 
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Here we present the first results on this phenomenon 
from a new 85 Rb BEC machine [18 j . finding good agree- 
ment between the measured collapse times and those pre- 
dicted by a GP model. Although we use the same atom, 
our experiment has several important differences from 
the original JILA work. Most notably, our condensates 
are confined in a purely optical potential, with a homo- 
geneous magnetic bias field applied to manipulate the 
interatomic interactions. In addition, we have measured 
condensate collapses with 4 x 10 4 atoms in a tighter trap, 
which together result in an initial density over an order of 
magnitude larger than in Ref. 0]. This leads to shorter 
collapse times and lower values of the critical scatter- 
ing length, but should not affect the ability of mean-field 
theory to describe the evolution of the system. It also 
allows us to investigate three-body recombination rates 
in a high density regime where they are the dominant 
source of atom loss. 

Our apparatus for producing Bose-Einstein conden- 
sates of 85 Rb with tunable interactions has been de- 
scribed in detail elsewhere [18] . In brief, we employ sym- 
pathetic cooling using 87 Rb as a refrigerant, initially in 
a quadrupole-Ioffe configuration magnetic trap and sub- 
sequently in a weak, large-volume crossed optical dipole 
trap. During the final evaporation, a magnetic bias field 
of 167 G is applied to reduce losses due to two-body in- 
elastic collisions |19l I20|. We can create condensates of 
up to 10 5 85 Rb \F = 2, m F = -2) atoms with a ther- 
mal fraction below 10% in a trap with harmonic oscilla- 
tion frequencies lo x ^ v , z = 2ir x {53, 22, 27} Hz. Conden- 
sates form at scattering lengths between a — +50ao and 
a = +200ao, where ao is the Bohr radius. We determine 
the scattering length from the applied magnetic bias field 
using the known parameters of the 155 G 85 Rb Feshbach 
resonance |21| . The field is calibrated by addressing ra- 
diofrequency transitions between the m f sublevels of the 
F = 2 manifold; the transition frequency is related to 
the magnetic field strength by the Breit-Rabi equation 
[2"2"] . The magnetic field can be determined in this way to 
within 5 mG, which near the zero crossing of the scatter- 
ing length corresponds to an uncertainty in a of ±0.2ao 

To observe condensate collapse, we follow the proce- 
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dure of Donley et al., tuning the atomic interactions using 
the Feshbach bias magnetic field as shown in Figure [TJ a). 
First, the scattering length is ramped smoothly from 
a = +89.1ao to a; n it = +0.2ao over 100 ms to produce a 
near-ideal, noninteracting gas. The magnetic field is then 
increased suddenly (< 100 fis) to a value at which the in- 



teractions are attractive a, 



collapse 



< 0, and held there for 



a time r before the trap is switched off and the scatter- 
ing length simultaneously increased to a = +50ao- The 
condensate is allowed to expand ballistically at this value 
for 15 ms, after which the magnetic field is switched off, 
changing the scattering length to czb g = — 443ao- Follow- 
ing a futher 5 ms of free evolution, the number of atoms 
present is determined by absorption imaging. 

The number of atoms remaining as a function of evo- 
lution time at a co n apse — — 20eto is shown in Figure [l[b). 
In agreement with the original work of Donley et aL, we 
observe a sudden and delayed onset of atom loss. This 
is explained by density-dependent three-body recombi- 
nation; when the interactions are made attractive, the 
condensate begins to contract slowly and its peak den- 
sity no increases, although not enough to cause significant 
three-body loss. As the condensate shrinks, however, the 
contraction accelerates, resulting eventually in a sudden 
implosion which increases the density by several orders of 
magnitude. This induces significant recombination losses 
(the loss rate scales with n 3 ) which ultimately halt the 
growth in density. The subsequent dynamics include fur- 
ther sporadic local implosions, which effect decay of the 
atom number in an approximately exponential form. We 
have observed remnant clouds surviving long after the 
collapse which contain several times the critical number 
of atoms N CT ~ 0.6ah o /|a|, where aho = \J h/muji w is the 
harmonic oscillator length [23 (the critical number for a 
condensate with a = — 20ao in our trap is N cr ~ 1200; 
c.f. Figure [T]). Such configurations have been shown to 
achieve stability through the formation of mutually re- 
pelling bright solitons [2]. 

The discrepancy between the JILA experiment and 
theoretical models concerns one of the most elemental 
characteristics of the bosenova: the time for which the 
atom number remains constant before the first density 
implosion - the so-called 'collapse time' i co iiapsc- Al- 
though the dynamics after the collapse are predicted to 
be complex and may exhibit behaviour beyond mean- 
field effects, the evolution prior to the first implosion 
should be captured in the mean-field approximation, and 
is determined almost exclusively by the initial density 
which is experimentally constrained. In particular, it has 
been noted that i C oiiapse does not depend strongly on the 
the three-body recombination rate K3 (SJ [TT] , which is 
not well-determined in the vicinity of the Feshbach res- 
onance. Despite this, Gross-Pit aevskii (GP) simulations 
were found to systematically overestimate the collapse 
time measured in the JILA system by almost 100%, a 
discrepancy at the 2a level given the experimental un- 
certainties |llj . 

As in previous work, we determine the collapse time by 
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FIG. 1: (a) Manipulation of the scattering length to induce 
and observe condensate collapse. After a variable evolution 
time r, the atoms are released from the optical trap simul- 
taneously with an increase in a from a co n apS e to +50ao. The 
cloud is allowed to expand at this value for 15 ms before the 
magnetic bias field is switched off. (b) Measured atom num- 
ber as a function of r for a C oiiapse = — 20ao- The solid line is a 
fit of the experimental data to equation |TJ . The atom num- 
ber remains approximately constant for a time t co n apsa , before 
a sudden onset of loss due to three-body recombination. 



fitting plots of the measured atom number versus time to 
the function 



N{t) = (N -N f )exp 
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+ N f , (1) 



for t > ^collapse, where No and Nf denote the atom num- 
ber at t < tcoiiapse and t > i co iiapse respectively. Figure [2] 
shows the collapse time determined in this way as a func- 
tion of a co iia P se for samples of No —Ax 10 4 atoms. As 
expected, the collapse time is shorter for larger |a co n apse |, 
as stronger attraction between the condensate atoms re- 
sults in more rapid contraction. The data are in qualita- 
tive agreement with the original experiment of Ref. [4 
and later theoretical work. 

To ascertain the ability of mean-field theory to quan- 
titatively reproduce our experimental data, we have per- 
formed numerical simulations for the parameters of our 
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system using the Gross-Pit aevskii equation for the con- 
densate wave function "F: 
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where Vt xap is the confining potential. Three-body recom- 
bination is modelled by the phenomenological inclusion 
of an imaginary loss term proportional to the three-body 
loss rate coefficient K$ (which differs by a Bose statistical 
factor of 3! from the coefficient for noncondensed atoms). 
This term leads to loss proportional to the cube of the 
atomic density: 
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and entails the assumption that the products of recombi- 
nation collisions leave the trap without interacting with 
the remaining atoms. As three-body processes dominate 
at the high densities relevant to this experiment [S] , we do 
not include the effect of two-body inelastic collisions. To 
make the computation tractable, in integrating equation 
([2]) we assume a cylindrically symmetric trap with oscil- 
lation frequencies u ZtP — 2ir x {53, 24} Hz, such that the 
mean trap frequency ui matches that of our crossed dipole 




FIG. 2: Collapse times as a function of scattering length for 
linit = +0.2a and No = 4 x 10 4 atoms. The data points rep- 
resent experimental values from measured decay curves such 
as that shown in Figure |TJb) , with error bars denoting the 
statistical uncertainty in the fit of equation Q. The solid 
line is the result of GP simulations for our experimental pa- 
rameters, and shows good quantitative agreement with the 
experimental data. The dotted lines represents the variation 
in the simulated collapse time due to experimental uncertain- 
ties in ai n it, a co n apso , u> and No. The dashed vertical line shows 
the critical scattering length for collapse at this atom num- 
ber, below which the condensate's kinetic energy stabilizes it 
against implosion. 



trap (the collapse time has been found to be relatively ro- 
bust to asymmetry in the trapping potential [11 ) . The 
simulation includes the 100 ms magnetic field ramp from 
a = +89.1ao to ai n it = +0.2ao, but we neglect the expan- 
sion of the condensate in our simulation, as the density 
spikes which trigger the recombination losses cease once 
the interactions are made repulsive. 

The results of this simulation for Nq = 4 x 10 4 and 
ainit = +0.2ao are overlaid with the experimental data 
in Figure [2] (solid line). For these simulations the three- 
body loss coefficient was scaled with a co n ap se as K 3 = 
8 x I0~ 14 a 2 cm 4 /s [13] • The dotted lines show the vari- 
ation in the simulated collapse time due to the combined 
experimental uncertainties in the initial scattering length 
and the trap frequencies, as well as run-to-run number 
fluctuations of 20%. The simulations show good quanti- 
tative agreement with the experimental data. 

It should be noted that the 100 ms ramp of the scat- 
tering length from a = +89.1ao to a = +0.2ao is not 
truly adiabatic. Our GP simulations show that the ramp 
excites breathing mode oscillations, despite the duration 
of the ramp exceeding the mean trap oscillation period 
by a factor of 3. The oscillation is predominantly along 
the weak trapping axes, and has an amplitude of ap- 
proximately 10% of the condensate radius. It occurs 
because, although a is varied smoothly, the condensate 
size does not depend linearly on the scattering length 
- in the Thomas- Fermi limit, the radius scales as a 1 / 5 . 
This excitation accelerates the contraction of the con- 
densate, decreasing t C oiiapse by approximately 15%. The 
effect is included in the simulations shown in Figure [2] 
It could be reduced by tailoring the magnetic field ramp 
to ensure that the condensate radius decreases smoothly. 
Our simulations show that ramping a 1 / 4 smoothly over 
100 ms reduces the breathing mode oscillations to below 
1%, and causes the collapse times to be indistinguishable 
from those for a condensate that is in the ground state 
immediately prior to the collapse. 

We now turn our attention to the possible systemat- 
ics which may affect the agreement between theory and 
experiment. The source of the largest experimental un- 
certainty in our system is the oscillation frequencies of 
the crossed dipole trap, which vary during the evapo- 
ration to BEC as the intensity of the trapping laser is 
reduced. For technical reasons, we cannot directly mea- 
sure the trap frequencies at the end of the evaporation. 
Instead, we make several measurements at higher inten- 
sities, which we fit to an analytic model of the dipole 
potential including the effect of gravity in the vertical 
direction. The model is further constrained by knowl- 
edge of the intensity Iq at which gravity overcomes the 
dipole potential and the trap vanishes. Due to the strong 
dependence of the vertical trap frequency on the laser in- 
tensity near Iq, and the variation in the intensity itself, 
we estimate an uncertainty in Q of 10%, predominantly 
in the vertical direction. As the peak density of a non- 
interacting condensate scales as w 3 / 4 , this corresponds 
to an estimated uncertainty in uq of 7%, which in turn 
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produces an uncertainty in the simulated collapse time 
of approximately 15%. This is not sufficient to explain 
the inconsistency between our results and the JILA ex- 
periment. 

We must also consider the possibility of a systematic 
error in our determination of atom number. No is cal- 
culated using the theoretical optical cross-section by in- 
tegrating the optical depth of an absorption image. We 
image the atoms on resonance with circularly polarized 
light and apply a small bias magnetic field along the 
imaging direction to provide a quantization axis. The 
calculation therefore makes use of the resonant cross- 
section and saturation intensity of the cycling transition 
\F = 3, m F = ±3) -> \F' = 4, m F < = ±4). As a result, 
our measured value of N Q is a lower bound: any errors 
in the polarization or detuning of the imaging light, or 
in the alignment of the quantization field, will reduce the 
measured atom number. We estimate the uncertainty in 
Nq due to these effects to be less than 5%. Furthermore, 
if the atom number were undercounted then correcting 
for this would decrease the simulated collapse times, as a 
higher initial density speeds up the contraction. This ef- 
fect therefore also cannot explain the disparity between 
our results and the original experiment, for which GP 
simulations overestimated the collapse times. 

Although £ co ii a p S e is only weakly dependent on the 
three-body loss coefficient K 3 , the shape of the loss 
curves is affected by varying this parameter. The values 
of K 3 used in simulations of the original JILA experiment 
ranged from K 3 = 2 x KT 28 cm 6 /s [S] to K 3 = 2 x 1CT 26 
cm 6 /s [IT]. Several authors also considered a relationship 
between the loss coefficient and the scattering length of 
the form K 3 ~ \a\ 2 for a < O [24], with Bao et al. 
deducing K 3 = 2.68 x lCT 13 a 2 cm 4 /s [13]. We find that 
these values cannot reproduce the shape of our measured 
loss curves. 

Figure [3] shows the results of GP simulations using val- 
ues of K 3 between 5 x 10~ 27 cm 6 /s and 5 x 10~ 29 cm 6 /s 
overlaid with experimental data for a co n apse = — 8.4ao. 
At higher loss rates, the high initial density of our sam- 
ple causes significant loss during the contraction of the 
condensate in the simulation during the moments leading 
up to the collapse. In fact, for K 3 > 10~ 27 cm 6 /s this 
initial loss is so great that there is no sudden implosion of 
the condensate and no discernible elbow in the loss curve 
|29| . In order to obtain the abrupt onset of loss that 
we observe in the experiment, a three-body loss rate of 
K 3 < 5 x 10 -29 cm 6 /s at a = — 8.4ao is required. From 
a similar analysis of the a co n apse = — 20ao data shown 
in Figure Qb) we find K 3 < 1 x 10~ 28 cm 6 /s at that 
scattering length. These limits are more than an order 
of magnitude below most of the values used to simulate 
the original experiment. Assuming a scaling with |et| 2 , 
they imply K 3 < 1 x 10~ 14 a 2 cm 4 /s. In this regime, loss 
after the initial implosion is caused by intermittent local 
density spikes between which three-body loss is negligi- 
ble. This was first predicted by Saito and Ueda jS] even 
before the JILA experiment. These discrete implosions 
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FIG. 3: Comparison of experimental and simulated collapse 
data for a co n apso = — 8.4ao. The data points show the mea- 
sured atom number N (normalized to No) as a function of 
evolution time r at a < 0, and the lines represent GP simula- 
tion results with different values of the three-body loss coef- 
ficient K3. A value of K3 < 5 x 10 -29 cm 6 /s is necessary to 
replicate the sudden onset of loss detected in the experiment. 



result in the numerous plateaus apparent in the simulated 
loss curve, although the scatter in our experimental data 
- caused primarily by run-to-run fluctuations in atom 
number - is too large to observe these directly. 

We have investigated the inelastic loss rates further 
by measuring the depletion of our condensates over time 
with positive scattering lengths, at which the condensates 
are stable. The rate at which atoms are lost due to two- 
and three-body inelastic collisions depends on the density 
profile of the condensate. In the limit that a — > 0, the 
density is given by the modulus squared of the ground 
state harmonic oscillator wavefunction, and the loss rate 
equation N/N = — J^i Ki{n i ~ 1 ) becomes: 

N = -N/t - r) 2 K 2 N 2 - rj 3 K 3 N 3 , (4) 

where r represents the one-body loss rate, r\ 2 = 
(2-Kal o y 3/2 and 773 = (v^Tra^)" 3 . In the Thomas-Fermi 
limit Na/dbn ^ 1, the condensate density takes on the 
shape of the confining potential and the loss rate equa- 
tion evaluates to: 

N = -N/t - -f 2 K 2 N 7 / 5 - j 3 K 3 N 9 / 5 , (5) 

with 72 = 15 2 / 5 /[147ra 3 / 5 a^ /5 ] and 73 
5 4 / 5 /[567r 2 3 1 / 5 a 6 / 5 a^ /5 ]. It should be noted that 
these expressions are valid only when the loss rate is 
small compared with the trap frequencies ^(n 1-1 ) <C w, 
so that the atomic density profile does not change 
significantly. 

Figure [4] shows the number of atoms remaining as a 
function of time in condensates with a = and a — 
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FIG. 4: Measurements of inelastic losses in 85 Rb condensates. 
The data points show the atom number as a function of hold 
time in the optical trap for condensates with a = and a = 
+37.6ao- The solid lines are fits of the solutions of p} and 
|5| to the experimental data, assuming K3 = (solid) and 
K2 ~ (dashed). Although the contributions of two- and 
three-body processes cannot be distinguished in this manner, 
these fits may be used to place upper bounds on the values of 
K2 and K3. 



+37.6ao. The lines plot the best-fit solutions to ([4} and 
([5]) respectively, assuming that the loss is entirely due to 
two-body (solid) or three-body (dashed) inelastic colli- 
sions. It is difficult to separate the contributions of two- 
and three-body loss purely from the shape of the decay 
curve, as has been noted in previous work Nonethe- 
less, attributing all of the measured loss to two- or three- 
body processes allows us to place an upper limit on the 
value of K 2 and K3 at these scattering lengths. Figure 
[5][a) shows these upper bounds for scattering lengths be- 
tween and +100ao. In our system, Na/a^ ~ a/ao and 
we use the Thomas-Fermi approximation ([5| except at 
a = 0. The error bars represent the statistical uncertain- 
ties in the fits; we assign an additional systematic error 
of 10% to incorporate the uncertainty in to. 

Theoretical calculations suggest that the recombina- 
tion rate should vary strongly with the two-body elastic 
scattering cross-section, with several authors predicting a 
universal K3 ~ a 4 scaling in the zero-temperature limit 
[2"5H2"T] . Our observations are consistent with a strong 
suppression of the recombination rate at the zero cross- 
ing of the s-wave scattering length, with the measured 
upper bound K 3 < (3.9 ± 0.7) x 10~ 29 cm 6 /s at a = 
an order of magnitude below that for a > +50a , and 
more than three orders of magnitude below the loss rate 
far from the Feshbach resonance, 2Y3 = 7 x 10~ 26 cm 6 /s 

We can combine this latest data with previous mea- 
surements of the two-body loss rate to further constrain 
the three-body recombination coefficient. Figure [5jh) 
shows the locus of possible K 2 , K 3 values for which the 
solution to Q best fits our experimental loss curve at 
a = 0. In Ref. [TU], Roberts et al. measured K 2 C — 
2.4 x 10" 14 cm 3 /s for thermal clouds in the vicinity of 
a = 0, corresponding to a value of K 2 — 1.2 x 10 
cm 3 /s for condensed atoms. This matches our measured 



upper bound of K 2 < (1.2±0.2) x 10" 14 cm 3 /s. Coupled 
with this result, our data is consistent with a value of the 
three-body loss coefficient K 3 < 10~ 30 cm 6 /s. Ref. [26_ 
predicts K3 ~ 5 x 10 -32 cm 6 /s at a = 0. In comparison, 
the three-body loss coefficient for 87 Rb is JC3 = 6 x 10~ 30 
cm 6 /s (28]. 

In conclusion, we have presented new experimental 
data on the collapse of 85 Rb Bose-Einstein condensates 
with attractive interactions in an optical dipole trap. Our 
results qualitatively match those of the original JILA 
bosenova experiment, but in addition agree quantita- 
tively with GP simulations. We find that a lower value of 
the three-body loss coefficient K3 than was used in sim- 
ulating the original experiment is needed to reproduce 
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FIG. 5: (a) Upper bounds on K2 (open circles) and Kz (filled 
circles) as a function of scattering length, calculated from fits 
to the solutions of Q and dBJ. The error bars represent sta- 
tistical uncertainties, (b) Locus of two- and three-body loss 
coefficients for which the solution to Q fits the experimental 
data for a — shown in Figure [4] The x and y intercepts 
correspond to the upper bounds shown in (a). Assuming 
Ki ~ 1.2 x 10~ 14 cm 3 /s [T5], the data suggest a three-body 



14 loss coefficient K 3 < lCT 130 cm b /s. 
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the sudden onset of loss that we observe. We have also 
analysed the decay of atoms from our condensates and 
thereby placed further constraints on the three-body loss 
coefficient K3 at small positive scattering lengths. We ex- 
pect that this work will inform future experimental and 



theoretical investigations of this rich quantum system. 
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